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ABSTRACT. This paper established that the mathematical modeling of Avian Influenza 
disease proposed has solution and unique when compare with other models in literature. 
Three equilibrium points were obtained, which stability of the system is investigated around. 
From analysis, the model stability depend on; first, 84 X must be less than (k+ n), secondly, 
the characteristic equation has all it coefficients of the A positive provided 

er+ted+rd> gS yO and Sy0(8» e + Bs a) 

Also, the model is stable at Avian Influenza Persistent equilibrium as the system satisfies 
J = F(0)G'(0) — F'(0)G(0) > 0 (Bellman and Cooke's [2]) condition. Hence, the model is 
locally asymptotically stable at both Disease Free and Persistent Equilibrium. 

Finally, the sensitivity analysis also, revealed that the factor responsible for the outbreak 
of Avian Influenza are majorly, four parameters and they are y, £5, £5, 0. These parame- 
ters are capable of preventing spread of Avian Influenza while restricting the movement of 
infected individuals during treatment session will go a long way manage the spread of Avian 


Influenza. 


1. INTRODUCTION 


Bird flu, also called avian influenza, is a viral infection that primarily infect birds, but 
humans and other animals are also susceptible, (World Health Organization- WHO [19]; 
MedlinePlus [13]; Bree and Sullivan [3]). It is rare for people to get infected with bird flu 
viruses, but it can happen. Two types, H5N1 and H7N9, have infected some people during 
outbreaks in Asia, Africa, the Pacific, the Middle East, and parts of Europe. There have also 


been some rare cases of other types of bird flu affecting people in the United States. Most of 
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the people who get bird flu have had close contact with infected birds or with surfaces that 
have been contaminated by the birds’ saliva, mucous, or droppings. It is also possible to get 
it by breathing in droplets or dust that contain the virus. Rarely, the virus has spread from 
one person to another. It may also be possible to catch bird flu by eating poultry birds or 
eggs that are not well cooked, MedlinePlus [13]. 

Bird flu illness in people can range from mild to severe. Often, the symptoms are similar 
to the seasonal flu, such as 

Fever Cough Sore throat Runny or stuffy nose Muscle or body aches Fatigue Headaches 
Eye redness (or conjunctivitis) Difficulty breathing 

In some cases, bird flu can cause serious complications and death. As with seasonal flu, 
some people are at higher risk for serious illness. They include pregnant women, people with 
weakened immune systems, and adults 65 and older. 

Treatment with antiviral medicines may make the illness less severe. They may also help 
prevent the flu in people who were exposed to it. There is currently no vaccine available 
to the public. The government does have a supply of a vaccine for one type of H5N1 bird 
flu virus and could distribute it if there was an outbreak that spread easily from person to 


person. 


2. LITERATURE REVIEW 


Mishra and Sinha [14] formulated mathematical model of Avian Influenza for both human 
and bird population. they computed basic reproduction number RẸ and RP for both hu- 
man and bird population respectively and locally and globally asymptotically stability were 
obtained for disease-free equilibrium and unique endemic equilibrium point was globally 
asymptotically state in bird population when R? > 1. In an extensive numerical simula- 
tions and sensitivity analysis for various parameters of the model were also carried out. the 
con concluded that critical analysis for the effect of Vaccination and Quarantined class with 
Recovered class. 

Putri, Watanabe and Nova [16] studied bird flu infection within a poultry farm with effect 
of spatial diffusion. Effect of threshold number that gives measure of infectious intensity 
was investigated, which depends on three factors; removal of infected bird, vaccination, and 
capacity of the farm. their analysis showed that transmission of bird flu within a poultry 


is expressible in terms of travelling wave solutions. It also shows a spatial significance of 


SENSITIVITY ANALYSIS FOR AVIAN INFLUENZA (BIRD FLU) EPIDEMIC MODEL 3 


diffusivity. These mean, infection can be fully controlled with appropriate vaccination and 
proper removal of infected individuals. 

Liu, Pang, Ruan and Zhang [10] discussed on global dynamics of Avian Infuenza epi- 
demic Models with psychological effect. The main consideration in the formulation was the 
psychological effect to Avian influenza in human population in which the Avian population 
exhibits saturated effect. The results demonstrated that the saturation effect within avian 
population and the psychological effect in human population cannot change the stability of 
equilibrium but can affect the number of infected human if the disease is prevalent. Also, 
the numerical simulation support the theoretical results and sensitivity analyses of the basic 
reproduction number gave effective control measures for avian influenze. 

Zhang, Li, Jin, and Zhu [21] presented an avian influenza A (HYN9) transmission model 
with vaccination and seasonality among human, birds, and poultry. the basic reproduction 
number for the prevalence of avian influenza is obtained. The global stability of the disease- 
free equilibrium and the existence of positive periodic solution were proven by the comparison 
theorem and the asymptotic autonomous system theorem. The study was finalized with 
numerical simulations to demonstrate the theoretical results. The simulation results revealed 
that the risk of H7N9 infection is higher in colder environment,and vaccinating poultry can 
significantly reduce human infection. 

Lee, Ko, and Jung [6] developed a spatial-temporal compartmental model that incorpo- 
rates the culling rate as a function of the reported farms and farm density in each town of 
the Republic of Korean. The epidemiological and geographical data of two species; chickens 
and ducks from the farms in the sixteen towns in Eumseong-gun and Jincheon-gun were 
used to find the best-fitted parameters of the metapopulation model. The local reproductive 
number RA using the next generation method is calculated as an indicator of virus trans- 
mission in a given area. Simulation results indicate that Ro has strong influence not only 
by epidemiological factors but also by geographical and demographical factors (density) and 
connectivity (or distance) farms. Based on the result obtained from simulation, they sug- 
gested the best culling radius of Preemptive (PE) culling should be adjusted by considering 


the local reproductive number in a targeted area. 


3. MODEL ANALYSES 


In view of the surveyed literature (Mishra and Sinha [14]; Putri et al., [16]; Liu et al., [7]; 


Zhan et al., [21]; Lee et al., [6]), it became imperative to investigate some properties of the 
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proposed models, specifically, equilibria states, threshold number for infection transmission 
and sensitivity analysis to know the most sensitive parameter possibly to cause bird flu epi- 
demics. Thus, subsequent subsections were devoted for the studies. The model is formulated 
based on the model proposed by Jabbari [7] and Liu [11] . The human population is classified 
into five subclasses; Susceptible individuals (S), Exposed individuals (E) Infected individual 
with avian strain (Ia), Infected individuals with mutant strain Im and Recovered individu- 
als (R) while the avian or bird population is classified into two subclasses; susceptible and 
infective denoted by X and Y respectively. 

The work of Jabbari [7] was extended by adding the exposed class and incorporating 
the progression rate of human from exposed to infections which is the parameter 0. Also, 
splitting infected human Jp used in Liu [11] into infected individuals with avian strain (Ia) 


and infected individuals with mutant strain (Im). 


3.1. Basic Assumption of Model. 
(1) Individuals are only recruited into susceptible sub population. 
(2) The number of susceptibles for the bird population is increased by new recruitment 
but reduced through natural death and infection. 
(3) The avian influenza virus is not contagious from an infective human to a subceptible 
human. It is only contagious from an infective avian to a susceptible human. 
(4) An infected avian keeps in the state of disease cannot recover, but an infected human 


can recover and the recovered human has permanent immunity. 


3.2. Model Equation. 
— =b- [AY se Bal, 22 B3Im)|S i uS 


"SERO [AY + Bola T BaIm] S "5 (0 + u)E 
dl, 





= 70E — (a + u + da + p1 Ma 


(1) —— =(1-y0E + al, — (qa + u + dm)Im 


dR 
dt 

dX 

H = c— b4 XY — kX 
Y 

as = ByXY — (k +n)Y 


= Yıla + Yəlm — uR 
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TABLE 1. PARAMETER DESCRIPTION 







































































PARAMETER | DESCRIPTION VALUE REFERENCE 

b Recruitment rate for human 30 Jabbari et al [7 

c Bird inflow 1000 Assumed 

£1 Rate at which bird to 0.2 Jabbari et al [7 
human avian influenza is contacted 

B» Rate at which human-human 0.4 Jabbari et al [7 
avian influenza is contracted 

b3 Rate at which human-human 0.0015 Jabbariet al [7] 
mutant influenza is contracted 

ba Rate at which bird 0.001 Jabbari et al [7] 
contract avian influenza 
Natural death rate for human 3.91 x 107? Liu et al [11] 

a Mutation rate 0.01 Gumel [5] 

da Deduced death rate due 1 Iwami et al. [6] 
to avian strain in human 

dm Induced death rate due 0.06 Jabbari et al. [8] 
to mutant strain in human 

k Natural death rate for birds 3.4246 x 1074 | Liu etal. [11] 

n Induced death rate due to avian strain in birds | 5 Iwami et al. [6] 

pı Recovery rate of human with avian strain 0.2669 Lucche et al. [12] 

p2 Recovery rate of human with mutant strain 0.05 Gumel [5] 

0 Disease progression rate of human 0.8 Assumed 
from exposed to infectious 

y Proportion of human 0.3 Assumed 
that are infected with avian strain 





3.3. Existence and Uniqueness of Model. Theorem 3.3.1 


(Derrick and Grossman, cited in Abdulrazak et al., [1]) Let D' denote the region |t — to| < 
a,||z — zo] € b,: £ = (zio, £20, .--, 240) and suppose that f(t,x) satisfies the Lipschitz 
condition ||f(t,z1) — f(t,z3)|| € k||zı — xe|| whenever the pairs (t,x ) and (t, x2) belong 


to R', where k is a positive constant, then, there is a constant 6 > 0 such that a unique 


continuous vector solution x(t) of the system in the interval t — to < ô. 


It is important to note that the condition is satisfied by requirement that 


be continuous and bounded in D’. 


We are interested in the region 


9f 


Ox,’ 





i=1,2,.. 
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0<a<R 











A bounded solution is sought ” D’ in the region whose partial derivatives satisfy 6 < a < 0, 
where o and 6 are positive constants. 

Proof: 

Let (1) be renamed as fı — f; 





fi) n —b- [BY + Bola + Balm|S — uS 
(fa) E = [EY + Bola + 851, ]S — (0+ u)E 
(fa) 2 = yE — (a + u+ da + 1)la 

(fj) dL (1 — 38E + ala — (p2 + H+ da) f 
(fs) a = (ils + qal, — uR 

(fe) “ = c — fAXY — kX 

(fr) am BaXY — (k +n)Y 
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Ofi| | Sofi Jof , 
E = | (Y + Bala + Balm + p) < 00; [D =0.< oo; [2f = | — 85S| < oo; 
Of fil | Jof =: Jof ; 
OL, = | — 83S| < oo; OR =e = 0< o; E =| f19| < oo; 
ð ð 
$h = [GLY + Bala + 63lm| < oo; | aeo ee | = |29| < os 
Of» Of» 
2E] = lus] < oes [ZB] = 0 < ool | = 0 < o0; [BË] = 1818] < oe 

(2) 
ô fo Ofe) _ fe) _ fsi ofo, _ 
ô fe fe 
I3x =|- ki colayl =1— BX] < oo; 
Ofr _ |f, Ov, oft) Ofr 
Pea aa = YS Oe laps sebs 
ð ð 
$E] = BY] < ooi $E] = [Bux — (k+n)| < oo; 


3.4. Region of Feasibility of Model Solution. Theorem 3.4.1 (For Host) 

The feasible region L defined by 

L= (co. E(t), Ialt), Int), R(t) €eR1:0€ N € t) is positively invariant and attracting 
with respect to the ODE for all t > 0 (Gilberto et al., [4]). 














Proof: 


To obtain the region in which the solution is bounded, the total population will be ; 


dN 7 ds dE di, dl, dR 
dt dt dt dt dt 
—b— pY S — Palas — BglmS — u$ + B1Y S + Polas + Ba1,,5 — 0E — wE+ 





«0 E — ala — ul, — dala — Yila + OF — YORE + oda — Gal — ula, — dala 


pila + Palm = LR 
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In the absence of induced death in the population (3) reduces to 





N 
OS psc | E | I, | TI. | R) 
dt 
«b— uN 
(4) 
Hence 














b 
0<N<-eER. 
H 


Theorem 3.4.2 (For Vector) 

The feasible region G defined by 

G= (axo, Y(t) eR2:0€ N «€ I is positively invariant and attracting with respect to 
the ODE for all t > 0 (Gilberto et al., [4]). 














Proof: 


To obtain the region in which the solution is bounded, the total population will be ; 


dN dX dY 
dt dt dt 
= c — pa XY — kX + pa XY — kY -nY 





In the absence of induced death in the vector population (3.4) reduces to 


dN 
—— «c—k(X -Y 
g k(X +Y) 
<c— kN 
(5) 
Hence 














0«N«zre R2 


3.5. Positivity of Solution. Theorem 3.5.1 

Given that the initial condition of the differential equations are Sp > 0, Eo > 0, Lag > 0, 
Imo > 0, Ro > 0, Xo > 0,Yo > 0 the solutions of S(t), E(t), Ialt), Im(t), R(t), X(t),Y (t) are 
non-negative for all t>0 (Gilberto et.al., [4]). 


Proof: 


Assume that 


t= sup{t > 0: S >0, E > 0, Ia > 0, Im > 0, R > 0, X > 0, Y > 0} € [0,0] 
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Hence t > 0; 


From the first equation in (1), 


(6) n > —(B1Y + Bola + B3lm)S > A > —(BÀiY + Bola + Bslm)at 


Integrating both sides of (6) gives 
dS 
(7) T >— i (81Y + Bola + B3Im)dt => In S(t) > —(6ıY + Bola + BaIn)t +K 


when £ = 0 in (7) implies 


(8) In S(0) = K 


Substituting back (8) into (7) yields 
In S(t) > —(BE1Y + Bola + B31,,)t + In S(0) 
=> In S(t) — In S(0) > —(f4Y + Bol, + B3lm)t 


Hence; 


(9) In —~ > —(81Y + Bola + B3lm)t 


PNY) ss p-(Y-BalatBsIm)t 


> S(t) > S (0)e- PY +BatatBatm)t 


But S(0) = So, hence 
(10) S(t) > Sye- Y iler Bsmt > Q 
Also from the second equation in (1), 


E E 
(11) > (0+ WB > SS > 0+ wat 
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Integrating both sides of (11) gives 


(12) [5 > - [+ ma ES eruta 


when t = 0 in (12) implies 


(13) In (0) = c 


Substituting back (13) into (12) yields 
In E(t) > —(0 + u)t + In E(0) 
=> In E(t) — In E(0) > —(0 + u)t 
Hence; 


E(t) 
(14) hm 2 6t 


Taking the exponent of both sides of (14) yields 


But E(0) = Eo,hence 


(15) E(t) > Ege OTHE > 0 


Considering the sixth equation in (1), 


dX 


(16) ps 


> (BY +h)X > > —(A4Y + k)dt 


Integrating both sides of (16) gives 


(17) a= / (GY + k)dt > In X (t) > —(E4Y 4- kyt +d 
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when t = 0 in (17) implies 


(18) In X(0) ^ d 


Substituting back (18) into (17) yields 
In X(t) > —(84Y + k)t + In X(0) 
=> In E(t) - In X(0) > —(84Y +k )t 


Hence; 


(19) In 





= X(t) > X(0)e-(&Y *t 


But X (0) = Xo,hence 


(20) X(t) ge Yt > 0 





Similarly, it can be shown that a(t) > I;ge-(^*^^*d«*e)* > 0, Im(t) > Impe2tetamt > 
0, R(t) UB UI Y(t) > he X-k- > 0 forall t> 0 


Hence, this completes the proof. 
3.6. Equilibria States. 


3.6.1. CASE I. In the absence of Bird flu and no infection i.e Y = 0, J, = 0 and Im = 0. 


The aim here is to get 
(21) e (Stu que FÉ R*, Xtym 


Q? Tn? 


To achieve (21), each compartments of equation (1) will be equated to zero in the absence 


of disease, that is; 
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b — [&1Y (t) + Bola(t) + Bstm(t)|S(t) — uS) = 0 
[BLY (t) + Bota(t) + Bslm(t)]S(t) — (u + 0) E(t) = 0 
YOE) — (a + u + da + q1)14(t) = 0 

(22) (1 — y)OE(t) + alalt) — (u + dm + v2) (t) = 0 
Prta(t) + Yalm(t) — R(t) = 0 
c— B4X (t)Y (t) — kX (t) =0 
B4X(t)Y (t) — (k -- n)Y (t) 20 


Now each variables in (22) will be solved as follows; 


From 6th equation 
(23) X= = 


and Human population dynamics is 
E = 0 from 2™ equation and R = 0 from 5” equation 


Hence, from 1st equation we have; 


(24) b— uS =0 


This implies 


b 
25 ue c 
(25) 7 
Hence 
b 
(26) (5 E, Ias Im R) - (2.0,0,0,0) } 
and 


m (em (60) 
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3.6.2. CASE II. Consider the situation whereby Y Æ 0, Im = 0, Ia = 0 


we use 














k+n 

28 X= 
(28) f, 
and from 6” equation 

C 
29 — 
(29) k + b4Y 
Equating (28) and (29) we obtain 

B4C — k(k 4- n) 

30 Y= 
0) Ba(k + n) 


Since J, = 0 and Im = 0, 1** equation reduces to 


(31) B-(u+fY)S =0 


This implies 


B 
32 Sr 
PA wt EY 
Substituting the values of Y into (32) yields 





(33) im (k 4- n)(uB4 — Bik) + 8184C 


and from 2" equation 
E 1 S 


(34) m 





After substituting Y, (34) yields 


T: B31(C84 — klk -- n)) 1 
(35) Bm To) rr? 





Thus, 


(36) Be = {9**, EDU 0f 


eg un Cpa- khan) 
(37) By = { ae } 





13 
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Remark:- It is obtainable to have bird flu in the vector population while human being only 
get exposed and this is inline with WHO, [20] report that it is hard to transmit the disease 


to human but when successful, mortality chance is 60%. Thus, we check for this possibility. 


3.6.3. CASE III. Persistence Equilibrium 
Suppose Y Z 0,1, #0 and Im 4 0 
Adding 1* and 2"4 equations yield 


(38) B — uS —(u+0)E =0 
Substituting (33) into (38) to get E yields 


—Bfhk(k +n) + Bi 84C B 
(ut+ 6)((k + m)(uBa + kB1) + 8184C) 


Expose leading to outbreak of Bird flu epidemic 


(39) E= 





From 3'7 equation, 





(40) ie Jn g —Bfik(k + n) + b1b4C B 
° a+pıtu+di (p+0)((k+n)(u8a kB) + 5C) 
(41) ga y0(bıBaCB — BBik(k + n)) 








(u + 8)(a + 1 + w+ da)(( + n)(uBa + kB) + B384C) 
Similarly from 4^" equation 
(1—»)0g9  ay@(GiBsaC — Byk(k --n))B — (ut dm + p2)Im 


42 = a = 
e) kz kı 1 P 





(1 — y)8gi(o + pı +u + da) — a70 (8184C — Pik(k+n))B 
kı(u + dm + 2) 





(43) E d= 
Obviously, from 5* equation, 
(44) R= seal; T1) 
where 

kı = (p. @)(at qi + u + da)((k + n)(uBa + kB1) + 8184C) 
(45) ko = (p + 0)((k + n)(uBa + kB1) + 6164C) 

gı = (- BBsk(k +n) + Bi 8aC) 
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3.7. Control parameter (Ro) for Model. The control parameter that governs the spread 
of disease is here obtained. The method of next generation matrix (Somma et al., {17]) is 
used. 


The model equation (1) is rewritten starting with the new diseases in the system of equation. 





“ = [BrY (t) + Bola(t) + Balm (0]8 (5) — (u + 0) E(t) 
2 = E(t) — (a + p+ da + 1) (0) 

(46) dl 
Gr = 0-7 VOE) + aL (5) — (u + da + 2) (5) 
as = B,X(t)Y (t) — (k + n)Y (t) 


From (46), f and v are deduced, where f refers to the new infections while v indicates other 


interactions in the infected compartments. This yields 





























[BY + Bola + Bala] (wt O)E 
B4 XY —y0E + (a + u + da + Y1)a 
T= and v = 
0 —(1 — y)8E — al, + (u + dm + ¢2)Im 
0 (k --n)Y 
Now differentiating f and v partially to yield F and V respectively, that is; 
Of: Of: Of əfi Ov Ov Av, Ou 
OE OI, Olm OY ðE la lm Y 
Of2 Of2 Of2 Ofe Ov? Ova Ovo Of 
_ | OE la lm ƏY _ | 3E la lm OY 
Pk Bhs A 5nd 9 = we Bin. e Bos 
OE Ol, Ol, OY OE la OL. OY 
Ofa Ofa Ofa Ofa Ov, Ov, Ov Ov, 
OE Ola lm OY OE la lm OY 
0 b25 B35 BS 
0 0 0 64S 
=> F= 
0 0 0 0 
0 0 0 0 
n0 0 0 0 
—y0 Bd, ocdqi 0 0 
and V= 
—(1—9)8 =a P2 + pt dm 0 
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Then 
1 
PET 0 0 0 
19 1 
(47) v-t = (u-4-0)a a 0 0 
0b a 1 p 
a(u+0)dd ad d 
0 0 0 E 
where 


a=pt+d,+a+y1 
b=yutydat+yy1-—a-p—da— yı 
and 


d = p2 + p + dm 


So the eigen value evaluated by |FV~! — AI| = 0 yields; 


12,34 — ( and A5 = i 


where 





g = Y h b2 — Y HU B3 + ^y Bad + ^y Bat — ^y Bada — Y Bat + et Ba + u B3 + Bsda + P391 
and 











j =p? +a uoo ud dro uaa 0d; --o80qps4- u34-u20-- 2d, uds, - uv uà pets dad 
LO dm+ pO pit pO potp dadm +H daqo +H dpi d- qo 0 dadm +0 dap2 +0 dnyit ip» 


where Ro = p(eigenvalues) 





Therefore, 


(48) Ry = —= 
The population will experience disease free if and only if Ro < 1 


This implies 


SO0g<j 


which implies that the force of infection must be less than the removal rate. 
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3.8. Stability Analysis for Model. In order to establish chances of controlling the dis- 
ease, the concept of Jacobian matrix will be used on the system. The Jacobian matrix of 


the model equation is obtained as; 
oF, OFF OF, OFF OF, OF, OF 
OS OE Ola Ol. OR OX ƏY 
OP, OP ôF, OF) OF) OF» OF? 
OS OE Ola Ol. OR OX ƏY 
OF; OF; OF; OF; OF, OF3; OFS 
ðS OE Ola Ol, OR OX ƏY 

J= OF, OF, OF, OF, OF, OF4 OF4 
OS OE Ola Ol, OR OX OY 
OF; OF; OF; OFs OFs OF& OFs 
OS OE Ola Ol, OR OX ƏY 
OFg OFe OF6 OF6 OFge OFe OF6 
ðS OE 0l, Ol, OR OX ƏY 
OF; OF; OF; OF; OF; OF; OFF 
0S OE Ola Ol, OR OX ƏY 


This yields 























—(m + p) 0 gus. -p35 0 0 -p18 

m —(0+u) B3S BS 0 0 fiS 

0 40 -d 0 0 0 0 
(49) J= 0 (1—4)0 a —e 0 0 0 

0 0 pı p2  —L 0 0 

0 0 0 0 0 “BY =- -64X 

0 0 0 0 BY  BX-(k+n) 
where 


m = p1Y + Bola T PBsIa 


d—octgituctd, e= p2 + H+ dm 





3.9. Stability of Disease Free Equilibrium for Model. At disease free, the 


Jacobian matrix (49) is considered in the absence of disease and hence reduce to 


=H 0 —B38 -3S 0 0 -1S 
m —(0+u) BS BS 0 0 fiS 
0 — «8 zc. a. vj. Dü 0 
(50) J= 0 (1—5)8 a —e 0 0 0 
0 0 pı p2 -u 0 0 
0 0 0 0^ "- —k -64X 
0 0 0 0 0 0 gX-(k+n) 
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Equation (50) is further reduced and the eigenvalue (J — AJ = 0) evaluation is carried out 


to obtain 
—(r+A) b2 5 p3 S As 
(51) TR y0 —(d + à) 0 0 255 
(1—3)8 a —(e+A) 0 
0 0 0 EX) =) 








Determining Stability at DFE of Model require evaluation of determinant of equation (51) 
such that the eigenvalues are negative or / and where necessary, characteristics equations 


results from the determinant satisfies Routh-Hurwitz Necessary Criteria, which stated that: 


(1) all the polynomial coefficients must have the same sign 


(2) all the polynomial coefficients must be non-zero 


In doing this, let q = 84X — (k + n) — A, and r = 0+ u, then equation (51) yields 
(—(r + A))(- (d + A)(-(e + 29) a — (828 y8(—(e) a) + 825804 —0 


Factorizing q gives 





a[-(r + A\(d+ (e +d) + 8306 + 8500] = 0 


This implies 





either 
(92) q=0 

or 

(53) —(r + A) (d A)(e - 3) + oS Oe + oS y0a=0 


Solving equation (52) yields 
(54) A= 84 X — (k 4- n) 


and solving (53) yields 


—erd—erA—edA—eM—rdA—rM -dX —M-FB,S40e-- Bo SyOA+ h3SyO0a=0 
This leads to 





Lm 
Or 
C 

Na 
Ll 

w 
> 
N 

—— 
c 
3 
a 


)+Aler+ed+rd—pb:9y0)+erd-—b:Sy0e-— bz Sy0a=0 
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From the solving, the system stability is concluded depending on; first, the equation (54), 
B4 X must be less than (k + n), and for the characteristic equation (55), which required 
examining the conditions of the roots. Thus, all the coefficients of the polynomial were 
found positive provided 
er+edt+rd> 85940 and Sy0(85e + 83a) 

Hence, the model is locally asymptotically stable at Disease Free Equilibrium. Next is the 


case of persistent of bird flu disease for the proposed model. 


3.10. Stability Analysis for Endemic Equilibrium for Model. In order 
to compute the stability analysis of the endemic equilibrium, equation (49) is considered in 


it full, but at simplification, it is reduced and the eigenvalue is evaluated as follows 


(56) 
—(m+oat+p)—-A 0 -b29 ps0 0 ES 
m —(0c-u)—^ oS B3 S 0 pı S 
" 0 79 ERSA TO 9 » 
0 (1-3)8 a —p-A " : 
0 0 0 0 | —AY —k-A —PaX 
0 0 0 0 BY BaX —(k+n)—2 
m =—(m+oa+A)(-(e+A))(-(4+ A))C- Ce + A) (p)+ 


(8 5 m 6(—(e + A)) (p) — (—8s 5 m y 6 a(p)) 


Factorising (p) 


(58) pi(b+ A)(c -- A)(d 4- A)(e - A) + Bo SmyO(e+ A) + B3 S my] 





This implies 


(59) p=0 


Or 





(60) (b+ A)(c 4- A)(d 4- A)(e +A) + £3 Sms 0(e -- À) - 63Smy0a=0 
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where 

m= pY F Bola + P3Im 
61) p — —(BAY +k t AY(B4X — (k +n) — A) + 82XY 
61 

c=p+0 

b=m+u 


Now evaluating (59) yields 





(62) 

21/9 V X28? —2XY Bj + Y?8? —2 Xnf, — 2Ynf4 +n? 
and 
(63) 





— 1/2 V X8? — 2 XY B + Y?6? — 2 XnB4 -2Ynf, +n? 


Also expanding (60) yields the characteristic polynomial 
(64) 
A +(b+c+d+e)à’+(bc+bd+be+cd+ce+de)à? + (bcd+bce+bde+cde 











+ BoSmy6)\ 4+ bcedet+ 85$ me f3Smy0o —0 


Using the concept of Bellman and Cooke’s Theorem [2] on equation (64) for establishment 
of the system stability. 
Theorem 3.10.1 
Consider the polynomial (64), which it is locally asymptotically stable around equilibria 
state if on the complex plane satisfies J = F(0)G'(0) — F'(0)G(0) > 0. 


proof: 

First, expressing the polynomial (64) in a complex plane by re-writing the polynomial con- 
taining A in terms if iy and Separate into real and imaginary as J(iy) = F(y) + iG(y) 
i.e. 

(65) 

J(iy) =(iy) + (b+ c-- d - e)(iy + (bc+bd+be+cd+ce+de)(iy)? +(bed+bce+bde 








+cde+f,oSmy6)(iy) - bcde-- 83S me-- £3 S m-y0o —0 
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J(iy) =y — (bc +bd+be+cd+cet+de)y?+bcde+ 83S me- B3Smy0a 





-i(( te+d+e)y*® + (bcd+bce+bdet+cde4 & Sm~6)y) =0 


Therefore, H (iy) = F (y) + iG(y), 
where, F(y) = y* — (be+bd+be+cd+ce+dejy +bcede+b:Sme+ b3Smyba 
and G(y) = ((b+c+d+e)y? + (bcd+bce+bde+cde4 B> S m»6)y) 
differentiating F(y) and G(y), then, substitute in the condition of the theorem, to have 














J = F(0)G'(0) > 0. Hence, the proof. Thus, the conclusion is that the system is locally 
asymptotically stable. 


4. SENSITIVITY ANALYSIS FOR MODEL 


Sensitivity analysis is carried out on each parameters, this is used in checking and identi- 
fying parameters that are liable to impacting the reproductive number. 
In order to obtain this technique,the concept adopted by (Tilahun, Makinde & Malonza, 


[18]) is here applied. The formula used in this concept for all basic parameters is defined as 


OR x 
Ro = Za. —* 
ey) A: ( Ox ) i Ro’ 


where x represents all the basic parameters. 


OR 
For instance, sensitivity index of Rọ with respect to y is AZo = ( =) x3 x — = 0.8683565952, 


Oy Ro 
all others are obtained in the same manner as Aj’? AUNT AD. AS, ATO ARo, AS. 


AR and are evaluated at y —0.3 0 = 0.8, pı = 0.2669, yo = 0.05, B2 = 0.4, B3 = 0.0015, 
u = 0.0000391, a = 0.01, 0 = 0.8, dm = 0.06, da = 1 The sensitivity indices is represented 


in table 2 
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TABLE 2. Sensitivity indices table 



































Parameter symbol | Sensitivity indices 
y 0.8683565952 
pı -0.1897545958 
(po -0.04201218574 
f» 0.9075403382 
b3 0.09245966216 
T -0.0001095245844 
-0.006800297646e 
0 0.00004887261139 
dm -0.05041462289 
da -0.7109576457 














4.1. Interpretation of Sensitivity Indices. In table 2, the parameters with positive in- 
dices i.e ( y, £5, £3, and 0 ) are those parameters that have great impact on the expansion of 
the disease in the community if their values are increasing. This is because the basic repro- 
duction number increases as their value increases, that is, the average number of secondary 
cases of infection increases in the community. 

Also from table 2, all the parameters in which their sensitivity indices are negative i.e (qi, 
P2, H, Q, dm and da) have capability of minimizing the disease in the community as their 
values increase while the others are left constant. As their values increase, the basic repro- 
duction number decreases, which leads to minimizing the endemicity of the disease in the 


community. 


5. CONCLUSION 


The sensitivity analysis revealed that factor responsible for the outbreak of Avian Influenza 
is majorly, four denoted parameters and they are y, £5, B3, 0. These parameters are capable 
of preventing spread of Avian Influenza while restricting the movement of infected individuals 


during treatment session will go a long way manage the spread of Avian Influenza. 
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